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ABSTRACT 

Jets and outflows from Young Stellar Objects (YSOs) are important signposts of currently 
ongoing star formation. In order to study these objects we are conducting an unbiased survey 
along the Galactic Plane in the 1-0 S(l) emission line of molecular hydrogen at 2.122/xm 
using the UK Infrared Telescope. In this paper we are focusing on a 33 square degree sized 
region in Serpens and Aquila (18° <1 <30°; -1.5°<b <+1.5°). 

We trace 131 jets and outflows from YSOs, which results in a 15 fold increase in the 
total number of known Molecular Hydrogen Outflows. Compared to this, the total integrated 
1-0 S(l) flux of all objects just about doubles, since the known objects occupy the bright end 
of the flux distribution. Our completeness limit is 3 • 10 _18 Wm -2 with 70% of the objects 
having fluxes of less than 10~ 17 Wm~ 2 . 

Generally, the flows are associated with Giant Molecular Cloud complexes and have a 
scale height of 25-30pc with respect to the Galactic Plane. We are able to assign potential 
source candidates to about half the objects. Typically, the flows are clustered in groups of 
3-5 objects, within a radius of 5pc. These groups are separated on average by about half 
a degree, and 2/3rd of the entire survey area is devoid of outflows. We find a large range of 
apparent outflow lengths from 4" to 130". If we assume a distance of 3 kpc, only 10 % of all 
outflows are of parsec scale. There is a 2.6 a over abundance of flow position angles roughly 
perpendicular to the Galactic Plane. 

Key words: ISM: jets and outflows; stars: formation; stars: winds, outflows; ISM: individual: 
Galactic Plane 



'■H 1 INTRODUCTION 



The interstellar medium (ISM) in galaxies is radically influenced 
by star formation. Giant Molecular Clouds (GMCs) are heated and 
excited by outflows from protostars and radiation from high-mass 
young stellar objects (YSO). Changes in chemistry and probably 
the turbulent motion in GMCs are a result of star formation, espe- 
cially massive star formation. Therefore, it is of great importance 
to understand the formation of stars. 

Jets and outflows from YSOs are sign-posts of currently ongo- 
ing star formation (e.g. Bally et al. (1995); Eisloffel (2000); Froe- 
brich & Scholz (2003); Davis et al. (2009)). Previous studies of star 
forming regions like the Orion A molecular ridge by Davis et al. 
(2009) and Stanke et al. (2002), DR21/W 75 by Davis et al. (2007), 
as well as the Taurus-Auriga-Perseus clouds by Davis et al. (2008), 
have shown a large number of outflows from YSOs. However, there 
are a number of open questions which are still to be addressed. For 
example: Is a large number of jets and outflows a common occur- 
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rence in other star forming regions (low and high mass)? Is star 
formation triggered in infrequent bursts or is it an ongoing, multi- 
ple epoch process in each GMC? The presence of jets from YSOs is 
an indication of a young population and active accretion while the 
sparsity of them in regions with a sizeable population of reddened 
sources shows a more evolved region with a larger population of 
pre-main-sequence stars. The dynamical age of a protostellar out- 
flow is 10 to 100-times less than the turbulent lifetime of a GMC. 
Therefore, the presence of a large number of outflows will be an in- 
dication of currently ongoing or multiple epochs of star formation. 

Outflows from YSOs are also a direct tracer of mass accre- 
tion and ejection and can be used to estimate star formation effi- 
ciency from region to region. This is particularly the case in high 
mass star forming regions where the efficiency is grossly affected 
by existing massive young stars, which influence the environment 
via their hugely energetic winds and intense UV fluxes. Further- 
more, in massive star forming regions, where young stars form in 
clusters and massive stars influence their lower-mass neighbours, 
photo-evaporation and ablation of protostellar disks can suppress 
accretion. As a result the mechanism that drives the jets and out- 
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flows is switched off. To what degree do these interactions affect 
accretion in YSOs? 

The position of individual protostars can be detected from jets 
and outflows while their evolutionary stage (e.g. Class 0, Class I) is 
closely related to the brightness in H2 emission (Caratti o Garatti 
et al. 2006). The mass infall/ejection history can be determined 
from jets since there is a correlation between the jet parameters, 
mass infall rates and accretion luminosities (Beck 2007; Antoni- 
ucci et al. 2008). 

Studies by Eisloffel et al. (1994) and Banerjee & Pudritz 
(2006) suggest that outflows are aligned parallel with the local 
magnetic field and perpendicular to the chains of cores. However, 
existing observations give mixed results without a definite answer 
(Anathpindika & Whitworth 2008; Davis et al. 2009). Is there a 
correlation between the mass of the driving source and the orienta- 
tion of the flow with respect to the cloud filament (i.e. are massive 
outflows more likely to be orthogonal to the filaments than flows 
from low mass stars)? Furthermore: What fraction of outflows and 
jets are collimated, and what fraction are parsec-scale in length? Is 
there a correlation between the median age of the embedded pop- 
ulation and the mean flow length? Are outflows sufficient in num- 
bers and energetic enough to account for the turbulent motions in 
GMCs? 

In order to answer the above outlined questions we need to 
have a representative sample of jets and outflows from young 
stars which is free from selection effects. This will allow us to 
perform a statistically meaningful investigation of the dynamical 
processes associated with (massive) star formation. Therefore we 
are conducting an unbiased search for jets and outflows in the 
Galactic Plane using the UKIRT Wide Field Infrared Survey for 
H 2 (UWISH2 - Froebrich et al. (2011), hereafter Fll). The data 
is taken in the H2 1-0 S(l) emission line at 2. 122/im, and thus 
highlights regions of shock excited molecular gas (T ~ 2000 K, 
nH 2 > 10 3 cm~ a ), as well as fluorecently excited material. Hence, 
it can be used to trace outflows and jets from embedded young stars. 

In this project we focus our attention on the Serpens/ Aquila 
region in the Galactic Plane, covered by UWISH2. In particular we 
investigate the area 18°<1 <30°; -1.5°<b <+1.5°, which approx- 
imately covers 33 square degrees (see Fig. 1 for an overview of the 
region). It is the first continuous area of this size completed and 
covers about 20 % of the total UWISH2 area. Hence, it will allow 
us to obtain first but statistically significant results about the popu- 
lation of jets and outflows from young stars in the Galactic Plane. 

In this paper we present our analysis of the above mentioned 
region. We focus on the detection of the jets and outflows, as well as 
their potential sources. We further analyse the spatial distribution of 
the discovered objects, their apparent lengths, position angles and 
fluxes. To convert apparent measurements (such as size) into phys- 
ically meaningful measurements we need the distance. We assume 
that all outflows are at a distance of 3kpc throughout the paper, 
since this turns out to be the most commonly measured distance 
for outflows in our sample. In our forthcoming paper (Ioannidis 
& Froebrich, in prep., Paper II) we will discuss in detail our dis- 
tance measurement method of the individual jets and outflows. We 
will then determine in detail statistically corrected luminosity func- 
tions and length distributions of the jets and outflows and investi- 
gate their distribution in the Galactic Plane. Finally, in Ioannidis & 
Froebrich (in prep., Paper III) we will investigate the driving source 
properties as well as the environment (cloud structure, clustered 
and isolated star forming regions) the jets and outflows are in. 

In Sect. 2 we discuss our data and analysis methods. We then 
present our results and discussion in Sect. 3, which includes the 



positions and distribution of outflows, the lengths and flux distri- 
butions as well as their orientation. We conclude our findings in 
Sect.4. 

2 DATA AND ANALYSIS 

2.1 Near infrared UKIRT WFCAM data 

We obtained our near infrared narrow band imaging data in the 
1-0 S(l) line of H 2 using the Wide Field Camera (WFCAM 
- Casali et al. (2007)) at the United Kingdom Infrared Tele- 
scope (UKIRT). The camera consists of four Rockwell Hawaii-II 
(HgCdTe 2048 x 2048) arrays with a pixel scale of 0.4". The 1- 
0S(1) filter is centred at 2.122 ^m with AA = 0.021 ^m. Our data 
is part of the UWISH2 survey (see Fll for details) and the H2 
images are taken with a per pixel integration time of 720 s under 
very good seeing conditions. The typical full width half maximum 
(fwhm) of the stellar point spread function is 0.7", the 5 a point 
source detection limit is about 18 mag (in broad-band K) and the 
surface brightness limit is about 10~ 19 Wm~ 2 arcsec -2 when av- 
eraged over the typical seeing (Fll). Our narrow band data were 
taken between 31st of July 2009 and 9th of September 2010. 

Furthermore, we utilised the UK Infrared Deep Sky Survey 
(UKIDSS) data in the near infrared JHK bands taken with the 
same telescope, same instrumental set-up and tiling as part of the 
Galactic Plane Survey (GPS, Lucas et al. (2008)) in order to per- 
form the continuum subtraction (H2-A') of our narrow band images 
and to generate three band colour images (JHK, JHH2, JKH2; 
see Sect.2.3). When compared to our narrow band data the NIR 
broad band data has very similar quality (fwhm and depth). Both 
data sets are taken typically about 2.5 to 4.0 years apart. 

Data reduction and photometry for both NIR data sets are done 
by the Cambridge Astronomical Survey Unit (CASU). Reduced 
images and photometry tables are available via the Wide Field As- 
tronomy Unit (WFAU). The basic data reduction procedures ap- 
plied are described in Dye et al. (2006). Calibration (photometric 
as well astrometric) is performed using 2MASS (Skrutskie et al. 
(2006)) and the details are described in Hodgkin et al. (2009). 

2.2 Difference Images 

In order to reveal mutually exclusive emission line regions asso- 
ciated with H2 jets and outflows (or other H2 emitters), the H2 
narrow band images have been continuum-subtracted using the K- 
band images from the GPS. Note that we are working with images 
about 13.3' x 13.3' in size throughout the project, since these are 
the images delivered by CASU. Every set of images (H2, K) is 
aligned based on their World Coordinate System (WCS). Photom- 
etry is performed to determine the mean fwhm of each image. 

To compensate for the different seeing conditions in both im- 
ages and to enhance the signal to noise ratio in the difference im- 
ages we Gaussian smoothed both images before the continuum sub- 
traction. The image with the worse seeing (image b) is smoothed 
using a Gaussian with a fwhm of r!j = 0.4". The image with the 
better seeing (image a) is smoothed with a Gaussian of a larger 
fwhm (rj), which we calculate as: 

r a s = yj{jwhm b y 2 - (fwhm a ) 2 + (ri) 2 

where fwhm a is the mean fwhm of image a and fwhm b is 
the mean fwhm of image b. Typically the seeing in both images is 




Figure 1. Positions of outflows (black squares) over plotted on a relative extinction map (grey scale image) based on near infrared colour excess in GPS data. 
Please see text for details on the background map. 



very similar and the smoothing radii and r b s are thus not very 
different. 

In order to completely remove continuum sources such as stars 
in the difference images, the if -band images need to be scaled. 
Since the scaling factor sc depends on the foreground extinction to 
the stars, we expect this to vary significantly across images taken in 
the Galactic Plane. We thus use a position dependent scaling factor 
sc(x, y) ('scaling image') for the continuum subtraction. 

For this purpose we determine the mean scale factor for all 
stars in small sub-regions (10" x 10"). Before applying the scaling 
image to the K-band it is smoothed with a radius of 30". Thus, 
we trace as much structure of the foreground extinction as possible 
while ensuring a sufficient signal to noise of the scale factor. 

The final difference images are then created as: 

Difference^ _x = H2 — K * sc 

Due to the Gaussian smoothing to adapt the stellar fwhm, the 
resulting difference images have a higher signal to noise than the 
original H2 and K images. While most of the non saturated stars 
are not present the H2 emission features are preserved. We note that 
the difference images are only used for the purpose of detection of 
jets and outflows from protostars (see Subsection 2.4) and not for 
photometry. 

2.3 NIR Colour images 

The search for H2 emission features (as described in Section 2.4) 
is mainly based on the K2-K difference images. However, the ver- 
ification that the observed feature is indeed a jet or outflow from 
a young stellar object is greatly facilitated by the additional use 
of near infrared colour composite images obtainable from our data 
sets. These images are also extremely useful to identify potential 
driving sources for the discovered jets and outflows. 

We thus created full resolution near infrared JHK, JHH2 
and JKR2 colour composites for every image. Each of the colour 
combinations enhances different aspects of the spectrum and there- 
fore makes it easier to visually detect and/or verify objects of in- 
terest. More specifically, the JKH2 images can be used to detect 
pure H2 emission regions. Since the H2 and K images have been 



taken at two different epochs, very 'red' or very 'green' objects can 
be identified as if-band variable stars in the JKH2 composites. 
In contrast, the JHK and JHR2 images can be used to detect ob- 
jects with if -band excess emission, in other words candidate young 
stellar objects. 

2.4 Outflow detection 

Our aim to investigate an unbiased sample of jets and outflows from 
young stars in the Galactic Plane cannot just be achieved by per- 
forming an unbiased survey. We also need to take great care not 
to introduce any detection bias when identifying the objects in our 
data. Thus, we follow the strategy described below to find all po- 
tential jet and outflow candidates. 

All H2-K difference images have been visually inspected for 
extended H2 emission features in full resolution. The order of in- 
spection was completely random to avoid the introduction of biases 
due to our search pattern. The original search has been performed 
by just one of us. However, the subsequent cleaning of this input 
catalogue (see below) has been done by two people. 

Every detected emission feature that has been identified in the 
H2-K difference images at first has been confirmed in the corre- 
sponding H2 image to avoid the inclusion of image artefacts. We 
then checked it against the full resolution JKK2 colour composite 
of the region (such as the one shown in Fig. 2). Thus, any object 
that could potentially be an image artefact has been removed from 
the source list. 

Finally, we need to clean the remaining objects from real emis- 
sion contaminants. These include e.g. Planetary Nebulae (PN), Su- 
pernova Remnants (SNR) and fluorescently excited regions such as 
edges of molecular clouds, HII regions and areas around young em- 
bedded clusters with massive stars. Planetary Nebulae can usually 
be identified by their visual appearance. In the case of SN shocks 
and fluorescently excited cloud edges this is more complicated, as 
they can mimic shocks from jets and outflows. We searched the 
SIMBAD database for known SNRs near to our identified objects. 
Any potential feature detected near such known SNRs has been ex- 
cluded from our list. We inspected regions with (obviously) fluores- 
cently excited cloud edges very carefully and excluded any object 
that potentially was not a jet/outflow feature. 
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The resulting list of H2 emission line objects is thus complete 
up to the survey detection limit, contains (almost) no false positives 
and is unbiased. 

2.5 Driving sources of outflows 

One vital task in order to understand the properties of the detected 
jets and outflows (e.g. the length distribution) is the identification 
of their potential driving sources. We utilise a number of published 
catalogues of YSOs as well as mid and far infrared source lists 
for this purpose. These include the catalogue of YSOs compiled by 
Robitaille et al. (2008) from the Spitzer GLIMPSE survey (Church- 
well et al. 2009), detections in the AKARI/IRC mid-infrared all- 
sky survey bright source catalogue (Ishihara et al. 2010; Yamamura 
et al. 2009), detections in the IRAS Point and Faint Source Cata- 
logue (Moshir 1989, 1991) and detections in the Bolocam Galactic 
Plane Survey (Aguirre et al. 201 1), which all cover our entire sur- 
vey field. Additionally we used Tf-band excess sources identified 
from JHK detected objects in the UKIDSS GPS, as well as K- 
band variable sources identified in the \i2-K difference images as 
potential driving source candidates. 

To decide on the most likely source for each H2 feature, we 
over-plot all potential source candidates over the K2-K difference 
images. Based on the vicinity of source candidates and H2 emis- 
sion, the alignment of a number of emission knots with a potential 
source, or the shape of bow-shock like features we grouped the 
H2 emission objects into outflows and assigned one MHO number 
to each of them, following the procedure outlined in Davis et al. 
(2010). Figure 2 shows MHO 2204 as an example of a detected jet 
with a potential source candidate. 

There are cases with several potential source candidates. If it 
is impossible to decide on one specific one we consider all pos- 
sibilities. In cases where non of the above mentioned catalogues 
allowed us to find a potential source for an H2 feature, we ad- 
ditionally searched the SIMBAD database for other indicators of 
YSO outflow sources, such as masers and (sub)-mm sources. Fi- 
nally, if there are several source candidates which apparently are 
the same object (such as an IRAS, Bolocam and Glimpse detection 
at roughly the same coordinates), we use as the source position and 
identifier the object in the survey with the highest spatial resolu- 
tion. There are a number of cases where we clearly detect a small, 
bipolar and symmetric jet but no object has been detected in any 
survey at the suspected source position. In these cases we assume 
the source to be situated between the H2 features on the flow axis, 
and it is given the identifier 'Noname' . Note, that MHO 2444 would 
be one of these objects, if it had not had deep JCMT data from Di 
Francesco et al. (2008). 

2.6 Photometry 

After the detection of the H2 emission features we have to measure 
their fluxes. This is straight forward, since our H2 images have been 
flux calibrated by CASU. 

We define apertures around each emission feature in the H2 
images. Great care is taken that continuum sources (such as fore or 
background stars at the same line of sight) are not included in the 
apertures. We also ensure that the apertures contain as little area 
as possible that seems to be free of emission in order not to add 
noise. Simultaneously, we use for each H2 feature a nearby sky- 
aperture that is completely emission free and defines the local sky 
level (see Fig. 2 for an example of the apertures used in the pho- 
tometry of MHO 2204). We then measure the number of counts in 



the H2 feature, corrected for the local background level and repeat 
that measurement, using identical aperture positions, in the scaled 
(using sc(x,y), see Sect. 2.2 above) A"-band image to correct for 
the continuum contained in the narrow band data. 

The background- and continuum-corrected H2 counts for 
each emission feature are then converted into physical flux units 
(W/m 2 ) using the flux zero points provided by CASU. The final 
flux values for the H2 emission in the 1-0 S(l) line (summed up 
over all knots in each MHO) are shown in Table A 1 in the Ap- 
pendix. 

Uncertainties in the photometry are based on the variation of 
the local background level in the H2 and if -band images and the 
choice of the exact position of the apertures enclosing the H2 emis- 
sion regions. Typical uncertainties of the measured fluxes are dis- 
cussed in Sect. 3.5. 



3 RESULTS AND DISCUSSION 

We have detected a total of 131 molecular hydrogen outflows from 
YSOs in our survey field of 33 square degrees. Out of these, 121 
(92 %) objects are new discoveries. Of the ten already known out- 
flows, two have been recently identified by Lee et al. (2011 in 
preparation) using UWISH2 data as well. Thus, 94 % of our out- 
flows are newly discovered in our data set. Hence, in the field in- 
vestigated here, the UWISH2 data increased the number of known 
molecular hydrogen outflows by a factor of 15. 

In Table A 1 in the Appendix we list the assigned MHO num- 
bers, positions, fluxes, apparent lengths, position angles, source 
candidates and their positions. In Table Bl we show K2-K images 
of each MHO and give a brief description of its morphology and 
potential sources. We note that there is a difference between the 
number of outflows (131) and the number of MHOs (134). The 
reason for this is that some of the already assigned MHO num- 
bers are part of the same outflow. More specifically MHO 2206, 
MHO 2207 and MHO 2208 are part of the same outflow. Similarly, 
the bow shock MHO 2212 is part of the same flow as MHO 2201. 
In cases like these we treat all MHOs as one outflow. All objects in 
Table A 1 that are part of an outflow containing several MHOs are 
marked with an asterisk. 

3.1 Spatial distribution and clustering of outflows 

We show the positions of all detected outflows in our field in Fig. 1. 
The grey-scale background map in this figure is a relative extinction 
map based on median near infrared colour excess (e.g. Rowles & 
Froebrich (2009)) determined from UKIDSS GPS data (Lucas et al. 
(2008)). Note that some regions in the map have small Ay off-sets. 
This is caused by known, minor photometric calibration issues of 
the GPS data which will be corrected in future data releases. Since 
we are not using the actual Av values in the paper, this is a pure 
'cosmetic' issue. 

As expected, the outflows are mainly located in areas with 
high extinction. Their overall distribution hence follows nicely the 
giant molecular cloud complexes visible in the Av map. Only a 
small number of objects is not situated within dense clouds. How- 
ever, there are numerous areas of high extinction where no outflows 
have been detected. In Paper III we will investigate in detail if there 
is an Av threshold for the detection of molecular hydrogen flows, 
and which fraction of high column density clouds does not show 
signs of on-going star formation in the form of jets and outflows. 

Histograms of the outflow positions along and across the 
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Figure 2. Left: Grey scale representation of the H2-X difference image of MHO 2204. The green ellipses indicate the apertures used to enclose the H2 flux, 
while the blue apertures represent the nearby sky. The solid line indicates the direction of the flow. Right: JKH2 color composite of the same region. Both 
images are 95" X 66" in size, North is to the top and East to the left. 



Galactic Plane are shown in Fig. 3. As one can also see in Fig. 1, 
there is a non homogeneous distribution along the Plane, with peaks 
at I = 20° and 1 = 25° and a minimum around I = 23° . This indicates 
the concentration of outflows to specific areas i.e. the giant molec- 
ular cloud complexes. The latitude distribution shows a Gaussian 
like distribution with a width of about one degree. Its centre is 
shifted to about b = -0.25° hinting that the main cloud complexes 
in this part of the Galactic Plane are at negative latitudes. If we 
assume a distance of 3 kpc for the outflows, the width of the dis- 
tribution corresponds to a scale height of about 25-30 pc. This is 
very similar to values for young stellar clusters (55 pc, (Friel 1995)) 
and OB stars (30-50pc, (Reed 2000; Elias et al. 2006)), but sig- 
nificantly smaller than the about 125 pc found for the dust at the 
solar galactocentric distance (e.g. Drimmel et al. (2003); Marshall 
et al. (2006)). The prominent increase of the number of objects at 
b = 1.5° (right panel of Fig. 3) is caused by a number of higher lati- 
tude clouds in particular at I = 26.5° (see also Fig. 1). 

We have investigated the clustering of the detected outflows by 
means of nearest neighbour distance (NND) distributions. Hence, 
we calculated the distance to the N th nearest neighbour for all out- 
flows in our sample and determined a histogram to analyse these 
distance distributions. The 1 st , 2 nd and 3 rd NNDs all show a sharp 
peak at small separations, indicating that the flows are typically 
found in close proximity to each other. This peak almost disappears 
in the 4 th NND distribution and is certainly gone in the 5 th NND 
distribution, where it is replaced by a wider peak at about 20-30' 
separation. 

Thus, we find that the outflows typically occur in groups with 
a few (3-5) members. The size of these groups is about 6' (a typical 
value for the 3 rd NND), which corresponds to about 5 pc if we as- 
sume all objects are at a distance of 3 kpc. This size is clearly larger 
than the typical embedded cluster (about 1 pc or slightly smaller; 
e.g. Lada & Lada (2003)) and smaller than typical nearby GMCs 
such as Taurus and Orion. Thus, currently ongoing star formation 
occurs in those 5 pc sized regions within GMCs and is not neces- 
sarily confined to embedded clusters. 

The groups of outflows are separated by about half a degree 
on the sky from each other. Given the number of members in these 



groups, the average group separation, the total number of objects 
and the field size, we can estimate that about two thirds of the sur- 
vey area are more or less devoid of molecular hydrogen emission 
line objects, a fact supported by the distribution shown in Fig. 1. 

3.2 Driving source candidates 

We have been able to assign driving source candidates, as described 
in Sect. 2.5, to 68 of our 131 molecular hydrogen outflows. This 
corresponds to just over half the objects. Overall, 75 source can- 
didates have been identified, since some outflows have more than 
one probable driving source. A complete list of the outflow driving 
source candidates is presented in Table Al in the Appendix. 

We find that typically brighter MHOs (according to their in- 
tegrated 1-0 S(l) flux, see Sect. 3.5 below) are more likely to have 
a source candidate assigned to them. In particular the average flux 
of MHOs with an assigned source candidate is almost five times 
higher than for MHOs without source candidate. For the median 
fluxes this ratio, however, is lower by a factor of three. In a future 
paper we will investigate if this is a pure selection effect. It could 
simply be that bright MHOs are close-by and hence their sources 
are easier to detect. On the other hand, bright MHOs could be in- 
trinsically luminous and thus be driven by brighter, easier to detect 
sources. 

3.3 Apparent outflow lengths 

For all outflows with assigned source candidate(s) we have mea- 
sured the apparent length(s). A histogram of the resulting lengths 
distribution is shown in Fig. 4. The light grey bars represent all 
the outflows (68), while the dark grey areas represent the already 
known objects (9). Note that the lengths are not corrected for the 
unknown inclination angles or the potential detection of just parts 
of the outflow. For single sided flows the lengths correspond to the 
separation of the source and the most distant H2 feature, while for 
bipolar flows the lengths correspond to the total length of the flow, 
a procedure also adopted in other surveys (Stanke et al. 2002). The 
known flows do not occupy a special part of the diagram, but rather 
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Figure 3. Histogram of the positions of the detected outflows along (left panel) and across (right panel) the Galactic Plane. 
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Figure 4. Histogram of the apparent outflow lengths. The light grey areas 
show all outflows, while the dark grey areas represent the known objects. 
In case there are several source candidates for one flow, the length corre- 
sponding to the most likely source is plotted. Note that an apparent length 
of 60" corresponds to about 0.9 pc for the assumed distance of 3 kpc. 



represent a random sub-sample. Note that in cases of several source 
candidates for an outflow, we plot the distance corresponding to the 
most likely source. The distribution will not change significantly if 
we plot any other lengths. 

There is a large range of apparent outflow lengths (from 4" to 
130") and a strong decrease of the number of outflows with increas- 
ing length. Of the 68 outflows, more than 60 % have an apparent 
length of less than 30" (or less than 0.4 pc if we assume a distance 
of 3 kpc). This is in agreement with the distribution of lengths found 
by Davis et al. (2009) and Stanke et al. (2002) along the Orion A 
molecular ridge and Davis et al. (2008) in Taurus-Auriga-Perseus 
(NGC 1333, L 1455, L 1448 and B 1). Only 10% of all outflows 
have a length of more than 1 pc (if they are at 3 kpc). In Paper II 
we will discuss the physical outflow length distribution in parsec 
after considering the individual distances to the objects. Note that 
the observed distribution does not agree with a randomly orientated 
sample of flows which all have the same length. 



3.4 Outflow position angles 

Previous studies (e.g. Mouschovias (1976)) suggested that clouds 
collapse along magnetic field lines to form elongated, clumpy fila- 
ments from which chains of protostars are born. Furthermore, the 
associated outflows are aligned parallel with the local magnetic 
field and perpendicular to the chains of cores (Banerjee & Pudritz 
2006). Our sample is ideal to study the relation between cores, fil- 
aments and outflows because of the large sample of flows that are 
distributed over a considerable number of molecular clouds. We 
hence determined the outflow position angles and list them in Ta- 
ble A 1 in the Appendix. 

A histogram of the position angles is shown in Fig. 5, using a 
bin-size of 30 degrees. The plot shows that a large fraction of out- 
flows has a homogeneous distribution of position angles, with five 
of the six bins occupied by 11.2 objects on average. However, at 
angles between 120° and 150° one can identify an over abundance 
of outflows. There are 20 outflows in this bin, which corresponds to 
a 2.6a deviation from the average of the other bins. We performed a 
Kolmogorov-Smirnov test to determine the probability that the ob- 
served position angle distribution is identical to a homogeneously 
distributed sample. For outflows with two potential sources, and 
hence two possible position angles, each angle was weighted by 
half; if there are three sources the weight for each angle was one 
third. We find a probability of just 10 % that our objects are drawn 
from a homogeneously distributed sample. 

We further investigated if there is a trend of objects with a par- 
ticular range of position angles with sky position. Nothing could be 
found. The objects falling into a particular bin in the position angle 
histogram are distributed completely homogeneously amongst the 
entire sample. 

The Galactic plane has an inclination of about 64 degrees with 
respect to the ecliptic in our survey area. Thus, the peak in the posi- 
tion angle histogram between 120° and 150° indicates that there is 
a 2.6a over abundance of outflows orientated almost (PA = 20° with 
respect to Galactic North) perpendicular to the Galactic Plane. It is 
therefore possible that our measured flow position angles are in- 
fluenced by the large scale cloud structure. Previous studies of the 
Orion A molecular cloud by Davis et al. (2009) and Stanke et al. 
(2002) and in Taurus-Auriga-Perseus by Davis et al. (2008) have 
shown a homogeneous distribution with no significant trends in the 
orientation of outflows. However, studies of DR21/W75 by Davis 
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Figure 5. Histogram of the distribution of position angles from the detected 
outflows in 30 degrees bins. The distribution shows a 2.6 sigma over abun- 
dance of flows with position angles between £ = 120° and I = 150°, or an 
orientation roughly perpendicular to the Galactic Plane. 



Figure 6. Histogram of integrated 1-0 S(l) fluxes (in logarithmic units) for 
the detected outflows. The light grey areas include all outflows while the 
dark grey areas represent the already known MHOs. Our completeness limit 
is about3-10 -18 Wm- 2 . 



et al. (2007) have shown that flows (in particular massive ones) are 
orthogonal to some degree to the molecular ridge. In Paper III we 
will analyse this distribution in more detail. We will measure, if 
possible, the orientation of the outflows with respect to its parental 
cloud filament, to verify or disprove the above findings. 

3.5 Outflow flux distribution 

The measured integrated 1 -0 S( 1 ) fluxes of the outflows are listed in 
Table A 1 in the Appendix. Figure 6 shows a histogram of this flux 
distribution. The light grey areas include all the outflows while the 
dark grey areas are the already known objects. Clearly, the known 
outflows dominate the bright end of the distribution and all but one 
are brighter than 10~ 17 Wm -2 . In particular, about 60% of the 
integrated flux in the 1-0 S(l) line of molecular hydrogen can be 
attributed to already known MHOs. Thus, our survey just about 
doubles the known integrated 1-0 S(l) flux from jets and outflows 
in the survey area, in contrast to the 15 fold increase in total outflow 
numbers. 

In general we find a sharp rise of the number of MHOs with 
decreasing integrated flux. More than 70 % of the MHOs have an 
integrated flux of less than 10~ 17 Wm -2 . The maximum number 
of MHOs occurs at about 3 ■ 10~ 18 Wm -2 which corresponds to 
about 10~ 3 solar luminosities in the 1-0 S(l) line of H2 at our as- 
sumed distance of 3 kpc. We hence conclude that this flux is the 
completeness limit for outflows detected in our survey. There are, 
however, a number of objects with a significantly smaller integrated 
flux, indicating that the actual detection limit for 1-0 S(l) fluxes is 
lower than even 10 _18 Wm~ 2 in some regions (most likely dark 
clouds devoid of fore/background stars). This is also supported by 
fact that many of the outflows with a small overall integrated flux 
consist of only one emission line feature, 

For each of the integrated fluxes we list the photometric uncer- 
tainties in Table Al. Typical uncertainties of the measured fluxes 
are of the order of 10%. For bright objects they are significantly 
smaller. Wherever the integrated flux is very low and the variation 
of the local background level is high we only determine an upper 
flux limit. 

How the measured integrated flux distribution converts into a 



statistically corrected luminosity distribution will be discussed in a 
forthcoming paper (Paper II). 



4 CONCLUSIONS 

In order to investigate the dynamical component of the star forma- 
tion process we perform an unbiased search for jets and outflows 
from YSOs along the Galactic Plane. Our data has been taken as 
part of the UWISH2 survey (Fll). It uses as a tracer the 1-0 S(l) 
emission line of H2, and here we focus our attention on a continu- 
ous 33 square degree sized region (18°<1 <30°; -1.5°<b <+1.5°) 
in Serpens and Aquila. 

We identify 131 outflows from YSOs from which 94% (123 
objects) are new discoveries in our data set. Therefore, our survey 
has increased the number of known molecular hydrogen outflows 
by a factor of 15 in the area investigated. 

We find a flux completeness limit for our outflow detection 
of 3 ■ 10~ 18 Wm -2 , with 70% of the objects showing fluxes of 
10~ 17 W m 2 or less. Typically, the already known outflows oc- 
cupy the bright end of the flux distribution. Our survey thus in- 
creases the known integrated 1-0 S(l) H2 flux from jets and out- 
flows only by a factor of two, compared to the large increase in the 
total number of flows. 

The overall spatial distribution of the detected outflows shows 
that they follow the GMC complexes in the Galactic Plane. Only a 
small number of objects is not situated within dense clouds. How- 
ever, there are large areas with high extinction where no outflows 
have been detected. Overall, about 2/3rd of the survey area is more 
or less completely devoid of jets and outflows. We further find that 
the flows typically occur in groups of 3 - 5 members with a size 
of about 5 pc (at our assumed distance of 3 kpc). These groups are 
typically separated by about half a degree on the sky. 

The distribution of flows perpendicular to the Galactic Plane 
shows a Gaussian like distribution with a width of about 25 - 30 pc 
(at our assumed distance of 3 kpc), similar to values for young stel- 
lar clusters and OB stars. 

We are able to assign possible driving sources to about 50 % 
of the outflows. Brighter MHOs are more likely to have a source 
candidate assigned to them. 
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We measure the apparent outflow length for outflows with as- 
signed driving sources. There is a wide range of lengths from 4" 
to 130" and a strong decrease of the number of flows with increas- 
ing length. More than 60 % of the outflows have an apparent length 
of less than 30" (or less than 0.4 pc if we assume a distance of 
3 kpc) while parsec-scale outflows are not common. Only 10 % of 
all outflows would have a length of greater than 1 pc at our assumed 
distance. 

The position angle distribution of flows with assigned source 
shows an 2.6 a over-abundance at angles between 120° and 150°. 
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APPENDIX B: MHO IMAGES 



Table B 1 : This table contains close-up grey scale YL2-K images for each detected 
MHO in our paper, as well as a brief description of the morphology of the emis- 
sion. We further briefly mention sub-structure of the objects as well as possible 
driving sources. All images are aligned with North to the top and East to the 
left. A scale is indicated in every case. In the top right corner we indicate the 
name of the H2 image which contains the object. The red line indicates the di- 
rection of the outflow while the red circles mark possible driving sources. Labels 
on each driving source candidate refer to the identification number in the source 
catalogue they are taken from (e.g. Glimpse, Bolocam, IRAS, GPS; see text for 
details). MHOs that are marked with a * after their MHO number are parts of one 
outflow. 
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